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Effect of anharmonicity of a cage potential for a magnetic ion vibrating in a metal is inves- 
tigated by the numerical renormalization group method. The cage potential is assumed to be 
one-dimensional and of the double-well type. In the absence of the Coulomb interaction, we 
find continuous crossover among the three limiting cases: Yu- Anderson-type Kondo regime, the 
double-well-type Kondo one, and the renormalized Fermi chain one. In the entire parameter 
space of the double-well potential, the ground state is described by a local Fermi liquid. In 
the Yu- Anderson-type Kondo regime, a quantum phase transition to the ground state with 
odd parity takes place passing through the two-channel Kondo fixed point when the Coulomb 
interaction increases. Therefore, the vibration of a magnetic ion in an oversized cage structure 
is a promising route to the two-channel Kondo effect. 

KEYWORDS: Kondo effect, impurity Anderson model, anharmonic potential, electron-vibration coupling, 
numerical renormalization group 



1. Introduction 

Recently, many researchers in condensed matter 
physics have been interested in characteristic ionic struc- 
tures, networks of cages filled or unfilled by guest ions. 
One typical example is the filled skutterudite compounds 
RT4Xi2(R = rare earth or alkaline earth; T — Fe, Ru, Pt 
or Os; X = P, As, Ge or Sb). In a certain case, the radius 
of filled ion is smaller than the size of cage. In such a 
case, it is expected that they will show various unusual 
physical behaviors due to the vibrations of the guest ions 
in strongly anharmonic potential. 

A very pecuhar feature is observed in SmOs4Sbi2.^"^^ 
It is reported that a large specific heat coefficient 7 is ob- 
tained, where the unusual phenomenon is its robustness 
against magnetic field. Some theoretical studies pro- 
pose as a possible scenario that local vibrations of the 
guest ions lead to the nonmagnetic Kondo effect. ^"^^ In 
these theories, the authors consider the situation where 
the guest ion moves back and forth among several po- 
tential minima in the cage potential. Actually, off-center 
modes of the guest ion are reported in several cage ma- 
terials.®"^^-* From a theoretical standpoint, there are two 
different types of local vibrations which couple with con- 
duction electrons; the breathing mode of the cage and 
the transverse mode, namely relative displacement be- 
tween the guest ion and the cage. We have studied the 
interplay between the electron correlation and the trans- 
verse type vibrations in Refs. 12 and 13. However, in the 
previous studies, the cage potential was assumed to be 
harmonic. 

Let us review former studies of the nonmagnetic 
Kondo effect concerning the transverse mode. Vladar and 
Zawadowski considered the model where a single atom 
tunnels between two positions, called as the two level sys- 
tem (TLS). They showed that the logarithmic divergence 
appears when the doubly degenerate states are connected 



each other through scattering processes of spinless con- 
duction electrons. ^''"^^^ This type of Kondo effect will be 
referred to as the double- well-type Kondo (DWK) in this 
paper. Subsequently, Yu and Anderson considered the 
first-order term of ionic displacement which produced the 
scattering processes between the spinless s-wave conduc- 
tion electrons to the p-wave ones."'^''-' When the electron- 
phonon coupling is strong, the ion displacement induces 
the doublet of polaron bound states even if the ion vi- 
brates in a harmonic potential, which may be called as 
the polaron doublet (PD). It was shown that the resul- 
tant effective potential for the ion displacement behaves 
like a double-well potential. Therefore, the model pro- 
posed by Yu and Anderson can be mapped to the TLS 
and is expected to show the nonmagnetic Kondo effect, 
which will be referred to as Yu- Anderson-type Kondo 
(YAK) effect. 

In this paper, we will study the effects of anharmonic- 
ity of a cage potential and discuss low-energy properties 
of a magnetic ion coupled with spinful conduction elec- 
trons. The cage potential is assumed to be of the double- 
well type, where the two minima are located symmet- 
rically around the center. By applying numerical renor- 
malization group method^^' to the present model, we 
find that in the noninteracting case, two types of non- 
magnetic Kondo effect mentioned above are realized and 
the low-energy properties make continuous crossover be- 
tween them when the shape of the cage potential is 
changed. Then, the role of Coulomb interaction U is 
investigated in the two distinct potential shapes where 
typical behaviors of the DWK and YAK effects are ob- 
served in the noninteracting case. We find that only in 
the typical YAK region, the 2-channel Kondo fixed point 
(2ch-K) appears with increasing U. This behavior is sim- 
ilar to the harmonic potential case. A notable point is 
that the YAK effect is more stable against the Coulomb 
interaction than the harmonic case. 
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2. Hamiltonian of the System 

When a magnetic ion with mass M vibrates in a cage 
potential V^age(Q)j the dynamics of the ion is determined 
from the following Hamiltonian, 

^^ion--|^(VQ)2 + yeagc(Q). (2.1) 

In the same way as the previous studies, we consider 
the situation that the ion vibrates in a one-dimensional 
potential. To investigate effects of anharmonicity of the 
cage potential, we expand Kage(Q) within the fourth- 
order of the ion displacement Q, 



(2.2) 



2^4 

where the cage potential is assiimed to be an even func- 
tion of Q and k2 and fc4(> 0) are the expansion param- 
eters. 

The shape of the cage potential is classified into two 
distinct types; single- well type for A;2 > and double- well 
(DW) type for k2 < 0. In this paper, we focus on the DW 
case. By introducing the dimensionless ion displacement 
q{= aQ), the Hamiltonian Hion is transformed into 
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Two potential minima are symmetrically located at ±qo 
in the DW potential and the height of potential barrier 
between them is Apot- 

The ionic eigenstates \rn) and eigenvalues are ob- 
tained from the Schrodinger equation of the ion system. 
We use the convention that m is numbered from zero 
in ascending order of the eigenvalues. Because of the re- 
Vcagei—Q), the eigenstates of the ion 

-q. The 



cage 
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|m) are characterized by the inversion, q 
operation of the inversion Vion is described by 
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Therefore, we define the parity of |m) as Pion = (1) 
depending on m = even (odd). 

We have already discussed the derivation of a general- 
ized impurity Anderson model for a magnetic ion vibrat- 
ing in a cage potential. Within the first-order of the 
ion displacement, the full Hamiltonian for a magnetic ion 
with an s-wave impurity electron orbital is written by 

H =Hc + Hhyh + ^^locai, (2.7) 
ffe =Y,e{k){cUk)co.{k) + cl{kMk)}, (2.8) 

H^yh =l^{K,cL(fc)/. + vAAk)Uq + h.c.}, (2.9) 

k(7 

/flocal =J^ion + £/ ^ fU. + Uflftflh, (2.10) 



where cocr(fc) (cicr(fc)) is an annihilation operator of the s- 
wave (j3-wave) conduction electrons and e{k) is assumed 
to be the dispersion relation with no angle dependence in 
the momentum space. For the localized impurity orbital, 
the annihilation (creation) operator is expressed by 
(/t) and e/ is its energy and U the Coulomb interaction. 
Note that the definition of Vi in the present paper is 
different from Refs. 12 and 13 by and identical to Vi 
in Ref. 13. 

3. Numerical Renormalization Group Approach 

To investigate the effects of the anharmonicity of the 
cage potential, we apply the numerical renormalization 
group (NRG) method^*' to the present model. The 
merit of the NRG algorithm is that it enables us to cal- 
culate the low-energy spectra and various physical quan- 
tities at finite temperatures with high accuracy in a con- 
trolled way. The key idea of the NRG is to discretize 
the continuous conduction bands in the logarithmic en- 
ergy scales characterized by A. The discretized Hamilto- 
nian is composed of the impurity site and the two Wilson 
chains. The two chains correspond to the ,s-wave and p- 
wave conduction bands, which are coupled with the im- 
purity site through the usual hybridization and the ion 
displacement- assisted one, respectively. The discretized 
Hamiltonian is written by 
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(3.1) 



where the n-th hopping matrix element is given by 

1 - 



^l_^-2n-l^l_^-2n-3' 



(3.2) 



and the matrix element qm,n and the Hubbard operator 
Xm,n are defined by 

gm,n = (m|g|n), (3.3) 
Xm,n = \m)®{n\. (3.4) 

Here, the density of states of the conduction bands, p, is 
assumed to be a constant, 1/2D, with the band width of 
2D. All the parameters with tilde are multiplied by the 
constant factor 2/{D(l + /l"^)}. 

We comment on the symmetries of the Hamiltonian 
(3.1). In the previous studies, three conserved quan- 
tum numbers arc used; the total electron number 7V*°*, 
the z component of total spin 5^°* and the total parity 
P. The last one is defined by the sum of the number of 
p-wave conduction electrons Np and the number of har- 
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monic phonons iVph, P = Np + A^ph = or 1 (mod 2). 
For the anharmonic potential case, obviously, we can use 
Pion instead of iVph- This straightforward replacement is 
based on the fact that the full Hamiltonian (2.7)-(2.10) 
is invariant under the inversion of the coordinate system. 

The Hamiltonian (3.1) is block-diagonalized and each 
block is characterized by the set of quantum numbers, 
A^*°*, Sl°^ and P. We treat the symmetric case, 2ey + [/ = 
0, and set various parameters for NRG calculations as 
follows; A = 3.0, band width D = 1.0 and M = 15000 
states kept at each NRG step. For the calculations of the 
ion eigenstates, there are two important parameters, the 
cutoff number for the states kept concerning the ionic os- 
cillations and that of the bases spanned by Hermite poly- 
nomials. We use sufficiently big cutoff numbers which de- 
pend on the potential shape controlled by go and Apot- 

As a result of the NRG calculations, we find that there 
are three types of the low-energy fixed points which are 
identified by analyzing the energy spectra. ^'^^ They are 
classified as the s-type, the 2-channel Kondo (2-chK) 
type and the p-type fixed points. Detailed discussions 
about the nature of these fixed points have been reported 
in Ref. 13. 

4. Crossover Behaviors among Three Different 
Regimes 

First, we discuss the results of the NRG calculations 
in the noninteracting case ([/ = 0) with Vb = 0.2 and 
w = 0.2 fixed. The middle panel in Fig. 1 shows the phase 
diagram in the parameter space of go and Apot ■ There are 
two curves represented by the circles and triangles. In the 
left part of the circle (triangle) line, no plateau of S'imp at 
fee log 2 is seen for Vi — 0.05 (0.10). The low-energy fixed 
point is always of the s-type in the entire phase diagram. 
Therefore, the two boundaries are crossover lines. 

Around the phase diagram, we show four graphs of the 
shape of the anharmonic potential with the low-energy 
levels represented by the dashed lines for (a) (go, Apot) = 
(0.5, 10) in the upper left, (b) (2.5, 10) in the upper 
right, (c) (0.5, 0.1) in the lower left, and (d) (2.5, 0.1) in 
the lower right, respectively. The diamonds in the phase 
diagram correspond to these sets of parameters. 

In Fig. 2(a), temperature dependence of S'imp is shown 
for various Apot with V\ ~ 0.05, C/ = and go = 2.0 
fixed. We may classify the parameter space into three 
characteristic regions depending on Apot by observing 
temperature dependence of Simp- In Fig. 2(b), Apot de- 
pendence of the thermal average of the square of the ion 
displacement (g^) at low temperatures is shown. 

We note first that for Apot < 1, the impurity entropy 
S'imp is released from the plateau at fee log 2 with the 
same temperature dependence as the usual Kondo effect. 
In this region, excitation energies of the ion do not show 
the character of the two local minima because the po- 
tential barrier between them, Apot, is low [see Figs. 1(c) 
and 1(d)]. In fact, the energy gap between and E^^^ 
is comparable to the difference between e'^^^ and e'^^^. 
Since the DW potential is shallow, the energy spectrum 
is not qualitatively different from the single well poten- 
tial. Figure 2(b) shows that the origin of the plateau is 
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Fig. 1. (Color online) Phase diagram in the parameter space of 
qo and Apot for Vb = 0.2, ZJ = 0.2 and U = (middle) and the 
shapes of the anharmonic potential with the low-energy levels 
represented by the dashed lines at (a) (go, Apot) = (0.5, 10), (b) 
(2.5, 10), (c) (0.5, 0.1), and (d) (2.5, 0.1). These four sets of pa- 
rameters are represented by the diamonds in the phase diagram. 
Note that the graphs (a)-(d) are drawn by different scales with 
regard to the q- and Apot-axes. In (b), the dashed lines with the 
symbol (2) represent almost degenerate energy levels. 

attributed not to the two potential minima induced by 
the anharmonicity of the potential but to the PD induced 
by combination of many excited states of the ion since 
the (g^) values are significantly enhanced over 4:{— Qq). 
Note that the distance from the center of the potential 
to one of the local minima is 2(= go). From these re- 
sults, it is reasonable to conclude that physical proper- 
ties of the impurity is described by the YAK effect. The 
role of the anharmonicity is to make the effect of the 
electron- vibration coupling stronger and the YAK effect 
is observed even for relatively small coupling constant 
Vi. 

When Apot is close to 1, no plateau appears in Simp- 
In this parameter region, the present system shows the 
characteristic behaviors that the s-wave Wilson chain in- 
cluding the /-orbital is perturbed by the p-wave Wilson 
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Fig. 2. (Color online) Temperature dependence of the impurity 
entropy Si^p for various Apot values with Vi = 0.05, U = 
and qo = 2.0 in the graph (a). The graph (b) shows the ther- 
mal average of the square of the ion displacement (ij^) at low 
temperatures. In the graph (c), Apot dependence of log]^Q(7/7o) 
is shown, where 7 is the specific heat coefficient calculated from 
the impurity entropy Si^p and 70 is that for Vi = U = case. 

chain through the electron-vibration couphng. From the 
analogy to the harmonic case, we call this physical situ- 
ation as a renormalized Fermi chain (RFC) regime. 

Lastly, when Apot is larger than 7, the energy gap 
between E^^^ and E^^l, AE'ion, becomes small, which is 
prominent in Fig. 1(b). Figure 2(a) shows that tempera- 
ture dependence of 5'imp for Apot > 7 is identical to that 
of the usual spin Kondo effect under a weak magnetic 
field. With increasing Apot, the origin of the plateau of 
5'imp at fee log 2 changes from the PD to the double- well 
potential minima. Correspondingly, Figure 2(b) shows 
that for Apot ^ 1, (9^) converges to 4. Therefore, the 
upper right region of the parameter space may be char- 
acterized by the double- well-type Kondo (DWK) regime. 

Actually, we can confirm realization of the DWK effect 
from the specific heat. Apot dependence of logiQ(7/7o) 
with go = 2.0 is shown in Fig. 2(c), where 7 is the specific 
heat coefficient and 70 is that for = [/ = case. In 
this study, 7 is calculated from central differences of the 
averaged impurity entropy S'imp('7i) at the temperature 
Ti corresponding to the z-th NRG step, 

j(J^'^ — 5'imp (ri) ~ '5'inip(ri+l) 

log^/A 

where 5imp(T'i) is defined by ^SimpiTi) + j{ShnpiTi-i) + 
5imp(r,;+i)). Ti is the averaged temperature on a loga- 
rithmic scale, v^Ti. The minimum of 7 is located in the 
RFC regime, where 7 is weakly enhanced by the electron- 
vibration coupling in comparison with the noninteracting 
Fermi liquid. The figure reveals that the picture of the lo- 
cal Fermi liquid is valid for the entire parameter region. 
The electron-vibration coupling leads to an effectively 
heavy local Fermi liquid state at low temperatures by 



the mechanism of the YAK for small Apot or the DWK 
for large Apot- 

Turning back to the phase diagram of Fig. 1, we com- 
ment on the Vi dependence of the boundaries determined 
from S'imp- For qq < 1.0 and Apot > 1, the boundaries 
determined from .S'imp have only a weak Vi dependence. 
This is consistent with the fact that the effect of the 
electron- vibration coupling is effectively weakened when 
the space where the ion can move becomes too small. 
Strong Vi dependence appears for qq > 1.0. The pa- 
rameter region where there is no plateau at fee log 2 is 
expanded with decreasing Vi. 

With regard to the double-well potential case, what 
type of the low-energy fixed point is realized is a subtle 
question. Let us start the discussion from the mapping 
of the present Hamiltonian (2.7)-(2.10) to the TLS. Ne- 
glecting higher ion excited states above the first one, we 
introduce the pseudo-spin to describe the almost doubly 
degenerate ionic eigenstates located at the two poten- 
tial minima. Concerning the ion state, the ground state 
with the even parity is chosen as the eigenstate of the z- 
component of the pseudo-spin, t]|0)= — |0). In the same 
way, r*|l)= |l). Similarly, we can define the pseudo-spin 
for the conduction electrons. For this purpose, we define 
oocr = fa and ai^ = J2k cicr(A:). 

Under this approximation, the operator of the ion dis- 
placement q corresponds to the x component of the Pauli 
matrix for the ion r^. Then, the Hamiltonian (2.7)-(2.10) 
can be rewritten as 

TT Tr—ST^ Sr^ t r ci i , ^^ion i 

a a,;3=0,l 

-f /fc + ^o5I{cLW«o.+h-C-}, (4.2) 

ka 

where r° is the x component of the Pauli matrix for 
the conduction electrons and g is the matrix element of 
q between |0) and |l). The second term represents the 
energy difference between |0) and |l), which originates 
from the tunneling process between the two potential 
minima. 

To begin with, we discuss the original TLS proposed 
by Vladar and Zawadowski'^'*"^^) for the spinless fermion 
case. By considering scattering processes of the second- 
order of the ion displacement, they showed that the TLS 
can be mapped to the anisotropic Kondo model under the 
pseudo-field given by Ai?ion- Because of the existence of 
the impurity /-orbital, the present model does not have 
the channel symmetry between the s-wave and p-wave 
conduction electrons^"^'* unlike the TLS. This asymmetry 
in the channel space is sufficient to realize the Kondo ef- 
fect for the present model even in the absence of processes 
of the second-order of displacement. Similar discussion 
is presented in the previous studies on the Yu- Anderson- 
type Kondo effect. ^^'^^^ 

Now, we go to the spinful case. Since the real spin can 
be used as the channel index for the conduction elec- 
trons, the non-Fermi liquid behaviors of the 2-chK may 
be expected for the spinful fermion case. It is known 
that the anisotropy of the (pseudo-)spin space is irrel- 
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Fig. 3. (Color online) Temperature dependence of the impurity 
entropy Simp for 0.0 < U < 0.4 in the upper left graph (a) and 
0.6 < U < 2.0 in the upper right one (b) with go = 2.0 and 
Apot = 0.08. That for Apot = 20 is shown in the lower graph 
(c). 



evant to the 2-chK fixed point. In the present NRG 
results, we have not seen any indication of the 2-channel 
Kondo effect for t7 = 0. We may argue that there are 
three reasons for the absence of the 2-chK fixed point. 
Firstly, AE'ion works as the transverse pseudo-field act- 
ing on the pseudo-spin. The non-Fermi liquid behaviors 
may be quenched by this external pseudo-field. Secondly, 
the eigenstates of the ion above the first excitation can 
not be neglected because the ion excitation energies are 
almost always lower than the band width D and plays 
the role of energy cutoff. Aleiner et al. showed that 
when higher ion excited states more than the first one 
are included into the TLS for the spinful fermions, the 2- 
chK temperature is always lower than Ai^ion- The Third 
reason is the channel asymmetry which is relevant to the 
2-chK fixed point. 

Lastly, we investigate effects of the Coulomb interac- 
tion U on physical properties of the impurity site. From 
the NRG calculations, we find that the responses against 
U depend on the potential shapes of the cage. Figure 3 
shows temperature dependence of the impurity entropy 
'S'imp for various U values with Vi = 0.05 and go — 2.0. 
In the upper two graphs (a) and (b), Apot = 0.08, and 
in the lower graph (c) Apot = 20. 

For Apot = 0.08, we observe the 2-channel Kondo 
fixed points a.t U = 0.201 and 0.709 from the analyses 
of the low-energy spectra. The type of the energy spec- 
tra changes from the s-type to the p-type at C/ = 0.201 
and return back to the s-type at 0.709 with increasing 
U. Such U dependence of the low-energy fixed points is 
the same as that for the harmonic case."'^^' In compar- 
ison with the harmonic case, there is one major quanti- 
tative difference about the lower critical U of the 2-chK 
fixed point. For the harmonic case, the lower critical U 
is roughly estimated to be the same order as the YAK 
temperature for the noninteracting case, see section 4.7 



in Ref. 13. This picture is not valid for the anharmonic 
case as shown in Fig. 3(a). It seems that the anharmonic- 
ity of the cage potential allows the YAK effect to be ro- 
bust against U, in particular for the shallow potential 
case. 

Completely different U dependence is observed for the 
case Apot = 20, as shown in Fig. 3(c). When U increases, 
a new plateau of S'imp at 2 fee log 2 appears. The value 
2fcB log 2 indicates that physical state of the impurity 
site is expressed by the direct product of the local mag- 
netic moment for the electron part and the two potential 
minima for the ion part. From the low-energy spectra, 
the fixed point is always of the s-type regardless of U. 
Temperature dependence of the entropy release from the 
plateau at 2 fee log 2 is controlled by the two energy scales 
corresponding to the s-channel Kondo temperature and 

Ai^ion- 

These two different types of responses against U cor- 
respond to the results of the noninteracting case. In the 
YAK regime, the physical state of the impurity is con- 
sidered to be the PD which consists of the spin singlet 
of the electron part and the linear combination of many 
excited states of the ion. By sharing the ion displace- 
ment by up and down spin electrons, there is an effec- 
tively attractive interaction, which leads to the energy 
lowering for the singlet channel. Introduction of repul- 
sive U changes the situation, and at the critical U, the 
singlet and triplet channels become degenerate. In this 
way, when U increases, the 2-chK fixed point is realized. 

On the other hand, in the DWK regime, the ion moves 
back and forth between the two potential minima. The 
ionic eigenstates are effectively written by only the al- 
most degenerate states and the higher ion excited states 
are not active. The real spin degrees of freedom define 
two independent screening channels for the pseudo-spin. 
However, in the noninteracting case, the 2-chK fixed 
point is suppressed by the three reasons mentioned pre- 
viously. When U increases, the up spin and the down 
spin electrons are correlated. Therefore, under the finite 
U, we can not use the real spin degrees of freedom as the 
independent channel index for the conduction electrons, 
and the 2-chK fixed point is not observed. 

5. Conclusions 

We have studied the generalized Anderson model con- 
structed for a magnetic ion vibrating in the double-well 
potential. By using the NRG method, we find that, for 
U — 0, the low-energy fixed point is always of the s- 
type and that the low-energy physics are described by 
the local Fermi liquid theory. The NRG results on S'imp 
and (g^) show that the parameter space can be classified 
into three different regimes: the YAK, DWK and RFC 
regimes. However, there is no sharp boundary between 
the three regimes and only smooth crossover behaviors 
are observed. 

One interesting question arises concerning the DWK 
regime. Since the effective Hamiltonian obtained by re- 
stricting to the two lowest vibrating states is represented 
by the TLS with the real spin as an additional channel in- 
dex, one may expect realization of the non-Fermi liquid 
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behaviors associated with the 2-chK fixed point. How- 
ever, the present NRG results do not show any indication 
of the 2-chK fixed point. We conclude that the renor- 
malization flow to the 2-chK fixed point is terminated 
by the following three factors: the transverse pseudo- 
field AE^ion, the higher ion excited states and the channel 
asymmetry inherent to the present model. 

Subsequently, we have considered the effect of 
Coulomb interaction U in the YAK and DWK regimes. 
In the YAK regime, by the same mechanism as the har- 
monic potential case, the 2-chK fixed point appears twice 
with increasing U. One interesting result is that the lower 
critical U of the 2-chK is considerably enhanced in com- 
parison with the harmonic case, even if the coupling con- 
stant Vi is not so large. We conclude that the vibration of 
a magnetic ion in an oversized cage structure provides us 
an opportunity to realize the 2-chK effect. On the other 
hand, the fixed point remains always of the s-type when 
U increases in the DWK regime. The present results indi- 
cate that the anharmonicity of the cage potential induces 
effectively strong electron-vibration coupling, especially 
in the shallow potential case. Therefore, wc may conclude 
that the anharmonicity of the cage potential is favorable 
to the realization of the nonmagnetic Kondo effect, which 
is one of the promising candidates for the magnetically 
robust heavy Fermion behavior observed in the skutteru- 
dite compound. 
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